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We discuss the sensitivity of a 100 TeV pp collider to heavy resonances produced in longitudinal 
vector boson scattering and decaying to a pair of Higgs bosons. A Monte Carlo study has been 
performed using the H ^ tt decay channel for both Higgs bosons, comparing the kinematics of such 
a signal to the irreducible Standard Model backgrounds. The results are presented in the context 
of a phenomenological model of a resonance {rj) coupling to goldstone modes, VlVl —>■ 9 —>■ HH, 
as can arise in composite Higgs models. With a fractional width of 70% (20%), the 5cr discovery 
reach is 4.2 (2.9) TeV in resonance mass for 10 ab“^ of integrated luminosity. We also discuss the 
dependence of the mass reach on the collider energy and integrated luminosity. 
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I. INTRODUCTION 

The discovery of the Higgs (H) boson [ 1 , 2] is one of 
the most exciting discoveries in physics. Its discovery 
completes the particle spectrum predicted by the Stan¬ 
dard Model (SM) [3, 4], and confirms the mechanism of 
spontaneous symmetry breaking to generate the longi¬ 
tudinal modes of the weak gauge bosons [5-11]. How¬ 
ever, the issue of very large radiative corrections to the 
mass of fundamental scalar particles has been acknowl¬ 
edged for some time [12]. This is the famous “natural¬ 
ness” or “fine-tuning” problem of the Higgs in the SM, 
where the quantum corrections due to quadratically di¬ 
verging loop integrals tend to drive the parameter values 
towards a very high energy scale. Beyond-SM theories 
containing additional symmetries (such as supersymme¬ 
try or other global symmetries) are able to protect the 
Higgs parameters from the quadratically-divergent radia¬ 
tive corrections. These symmetries imply the existence of 
new particles in loops which cancel or mitigate the loop 
corrections due to SM particles. 

In a class of models [13-15], the SM Higgs boson is 
itself another member of a set of Goldstone modes, gen¬ 
erated from the spontaneous breaking of a larger global 
symmetry. Other Goldstone modes are the and the 
Zp, which become the longitudinal components of the 
weak gauge bosons after electroweak symmetry break¬ 
ing (EWSB). In a subset of these models [16-18], a 
common feature is the spontaneous breaking of a global 
SO{f>) —>■ S'0(4) symmetry, where the four Goldstone 
modes are the four real components of the SU{2)l dou¬ 
blet Higgs field. 


The spontaneous breaking of the larger global symme¬ 
try is postulated to be caused by the formation of a con¬ 
densate, due to new strong dynamics at a high energy 
scale. As Goldstone bosons are automatically massless 
due to the Goldstone theorem, the question changes to 
why the Higgs boson mass is non-zero but small relative 
to the higher-energy compositeness scale. Various mod¬ 
els have been proposed to explain this “little hierarchy” 
between the EWSB scale v and the new compositeness 
scale fr). 

A signature of this strong dynamics in the Higgs sector 
would be resonances coupling to the Goldstone modes, 
i.e. the longitudinal W and Z bosons and the Higgs bo¬ 
son. The lightest of such resonances would preferentially 
decay to the lightest particles interacting with this sec¬ 
tor, ie. Wl, Zl, and H. As a benchmark, we study the 
phenomenological model for a new scalar resonance ry [19] 
whose Lagrangian is 

C = Csu + \d^-qd^r] - \mlrf' + ^r7i9'"7r“a^7r“ (1) 

^ ^ Jr) 

where 7 r“ represents the quartet of Goldstone modes in 
the Higgs doublet, plays the same role as the “pion 
decay constant” (/^) in QGD chiral perturbation theory 
and Qri is a dimensionless coupling. 

In composite Higgs models, rather than fully unita- 
rizing longitudinal gauge boson scattering, the Higgs bo¬ 
son delays perturbative unitarity violation until the scale 
frj. At this scale the strong resonances appear which to¬ 
gether unitarize the amplitude. In the scalar resonance 
model we look at, there is one resonance 77 below the com¬ 
positeness scale which partially unitarizes the amplitude. 
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depending on the coefficient a^. When = 1 the r] com¬ 
pletely unitarizes the amplitude, hence in this work we 
take = 1 to simplify the high-energy behavior while 
characterizing new strong dynamics in vector boson scat¬ 
tering (VBS) [20, 21]. Assuming no other decays other 
than to goldstones, in the high mass limit the width of 
the rj resonance is given by [19] 


r,= 


2 S 

Sirp 


( 2 ) 


implying that for our choice of a,; = 1 the fractional 
width is given by Trj/mjj = m^/(87r/^). 



FIG. 1. A display of a simulated event generated using the 
process pp —>■ r^jj with the decay chain p —>■ HH —>■ 4r at a 
100 TeV pp collider, with m,, = 3 TeV. Two Higgs-jets are 
shown with the large blue cones, each containing a pair of r 
leptons. Two small blue cones show the forward jets. The jets 
are reconstructed using the anti-fcr algorithm [22] using the 
Fast Jet package [23]. Yellow lines show charged hadrons. 


kinematic distributions associated with t] production and 
decay. The combination of the information garnered from 
measuring the different branching ratios and the associ¬ 
ated kinematic distributions can provide the definitive 
test of the Goldstone nature of the Higgs doublet field 
and information on its coupling to the new strong dy¬ 
namics. Sensitivity to double-Higgs production in the 
context of the SM and other theoretical approaches has 
been investigated [31-49]. The decay channels for the 
Higgs boson pair are indicated in Table VI, ranked by 
branching ratio [50] . In this table, we require that gauge 
bosons decay to leptons in order to suppress enormous 
backgrounds from QCD jet production. The first three 
channels in the table, HH —>■ 46, HH —>■ 262r, and 
HH —>• ivii'bb are subject to large QCD backgrounds 
from 6-jets and tt production respectively. These consid¬ 
erations motivate the studies using the 4r and 2^2b final 
states. The 2^2b final state has been studied both at 14 
TeV [51-53] and 100 TeV [39, 40, 54, 55] due to its good 
mass resolution. In this work, we present the first study 
of the 4t final state. The 2^2b final state may be further 
studied in the future as a complementary channel with 
different backgrounds, providing additional discovery po¬ 
tential. 


II. THE 4r FINAL STATE 


This analysis uses all decay channels for the r lep¬ 
tons from the Higgs bosons. An example of the decay 
T] ^ HH —>■ 4t is shown in Fig. 1. It shows two Higgs 
bosons with transverse momenta above I TeV, arising 
from the decay of the rj boson with a mass of 3 TeV. The 
decay of each Higgs boson leads to a single jet containing 
two T leptons. The event display was created using the 
Delphes fast simulation [56] and the Snowmass detector 
setup [57]. The event, generated with the Madgraph5 
Monte Carlo generator and showered with PythiaS [58] , 
was taken from the HepSim repository [59]. A complete 
color version of this event is shown in the Appendix. 


The sensitivity to the scattering process WpWp —>■ 
r] —>• WlWl at a 100 TeV pp collider has been stud¬ 
ied [24-30]. The W'^W~ channel faces substantial back¬ 
ground from tt production. It is interesting to probe p 
production using both the VlVl and HH final states. 
The model predicts branching ratios in the proportions 
2:1:1 for the WpWp, ZpZp and HH branching ratios 
based on counting the number of Coldstone modes in the 
Higgs doublet field. The distinctive feature is that Cold- 
stone modes have purely derivative coupling, which pre¬ 
serves the shift symmetry of the Coldstone fields. Specif¬ 
ically, the coupling terms in Eqn. 1 expand to pW^ , 
pZ'^Zf^L, and pd^Hdf^H respectively. The complete in¬ 
teraction term between the p resonance and the gauge 
bosons has the form p[2g'^Wl^W^^~ + (gW^^ — 

The specific Lorentz structure of these terms dictates the 


This study is performed without simulation of the de¬ 
tector response, but we include the effect of the r-tagging 
efficiency which is the most important performance char¬ 
acteristic. Hence it is likely that the inclusion of a realis¬ 
tic detector simulation with a similar r-tagging efficiency 
will not change significantly the results of this analysis. 
We assume that the requirement of four r leptons will 
suppress mis-identification backgrounds to the level that 
SM processes producing four prompt r leptons will domi¬ 
nate the backgrounds. This assumption is justified based 
on the r-lepton identification efficiency and QCD jet re¬ 
jection achieved by the LHC experiments. For instance, 
the hadronic decays of the r-lepton are identified with 
an efficiency of 60%, and with a QCD jet efficiency of I- 
2%, as reported by the ATLAS experiment [60-63]. The 
analysis most similar to our analysis is the high-mass 
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FIG. 2. Example of Feynman diagram for the process pp —>■ 
Vjj HHjj —>■ 4:T + jj via the production of the rj resonance 
in longitudinal vector boson fusion. 


Z°' —rr search, where the transverse momenta (px) of 
the r leptons are similar to our signal kinematics. In the 
double-hadronic mode, the dominant background arises 
from the 7*/^ — >■ tt Drell-Yan (DY) process, followed 
by multijet and VF/Z-|-jet background. The latter back¬ 
grounds are a factor of 3-4 smaller than the irreducible 
DY background. In the leptonic-l-hadronic decay modes, 
the DY, IF-l-jets, tt, diboson and single-top backgrounds 
all contribute approximately equally. In order to esti¬ 
mate the fake backgrounds and compare them to the ZZ 
background for our HH search, we consider the dibo¬ 
son analysis for ZV —>■ lljj [64]. This analysis shows 
that requiring Z ^ ll and V —^ jj with a mass win¬ 
dow cut yields a Z-bjets rate which is about 20-50 times 
larger than the VV rate. Given that the hadronic r se¬ 
lection is about 15 times more efficient for prompt t’s 
compared to QCD jets, the requirement of two additional 
t’s will suppress the Z-bjets background to a fraction of 
the ZZ background. Multi-jet background will be re¬ 
duced to a negligible level. Reference [64] also shows 
that the tt background is negligible in the high pt{V) 
region. The requirement of two additional r’s will sup¬ 
press diboson-bdijet and single-top-bdijet backgrounds to 
a negligible level. 

Dedicated studies with full simulation will be needed 
to design the future detectors which can maintain the 
T identification performance at high pt(t) at the same 
level that the LHC experiments have demonstrated. The 
above discussion shows that, if this performance can be 
achieved, the sensitivity studies presented below using 
the irreducible backgrounds should provide a reliable es¬ 
timate of the discovery potential. 


III. MONTE CARLO SIMULATIONS 


The goal of this analysis is to estimate the discovery 
potential of a future 100 TeV scale pp collider, based on 
the expected event rates and distributions for the signal 
and backgrounds after kinematic and fiducial cuts. The 


T-tagging efficiency, which is the main characteristic of 
detector performance, is assumed to be 60% similar to 
the LHC experiments [60-63]. 



mn (TeV) 



itin (TeV) 


FIG. 3. Cross sections for the vector boson fusion process 
pp —>■ pjj with p ^ HH —>■ 4r for the fractional resonance 
width of 20% (top) and 70% (bottom), computed using the 
MADGRAPH5 program at LO QCD. The following generator- 
level cuts have been applied: rrijj > 1.5 TeV and pT(jet) > 
50 GeV. 

The analysis was performed using the PythiaS [58] 
and Madgraph5 [65] MC models with the default pa¬ 
rameter settings. The MSTW20081o68cl [66] parton den¬ 
sity function (PDF) set was used. 

The signal cross sections as functions of are shown 
in Fig. 3. These cross sections scale approximately as 
powers of the resonance mass, m““, where 3.2 < a < 
4.8 (3.9 < a < 5.2) for a fractional resonance width of 
20% (70%) and i/s = 100 TeV. The range of a shows 
the departure from a constant-power law as 1 < < 

7 TeV, with the larger values of a corresponding to the 
larger values of m,;. Thus, the value of a can be used 
to estimate the slope of the curves in Fig. 3 at different 
values of m^. The dependence of the cross section on 
collider energy may also be parameterized as a power law, 
(•yi)^, with 1.8 < b < 3.3 (2.2 < b < 3.6) for a fractional 
resonance width of 20% (70%). Again, the range of b 
corresponds to 1 < nirj < 7 TeV, with the larger values 
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of b corresponding to the larger values of 

The HepSim public repository [59] was used to store 
simulated events in the ProMC file format [67, 68]. The 
samples were analyzed with a C++/ROOT program [69]. 
The jets were reconstructed with the anti-Zcr algorithm 
[22] using the Fast Jet package [23]. The typical choice 
of the distance parameter R for jet reconstruction is 
R ~ 0.4 for light-quark and gluon jets, and this value 
is motivated for the reconstruction of the forward jets. 
A smaller value of i? ^ 0.05 is motivated for the re¬ 
construction of the highly-boosted r-jets present in our 
samples. For simplicity, a single value of i? = 0.2 is used 
in this study, since the kinematic distributions are fairly 
broad and the discrimination between signal and back¬ 
ground is not sensitive to optimization of the choice of 
R. We assume that the development of sophisticated r- 
reconstruction algorithms and the use of jet sub-structure 
information in the ultimate data analysis will permit the 
separation of boosted —>■ rr and Z —>■ rr jets from 
mis-identified QCD jets and electrons mimicking r-jets. 

Jets have a generator-level requirement of prOet) > 
50 GeV, based on studies performed for the high- 
luminosity (HL-) LHC [70, 71] where this requirement 
was applied to suppress pileup jets in the forward region 
to an acceptable level. Given the hard dijet-mass (mjj) 
spectrum arising from the VBS topology, a generator- 
level requirement of rrijj >1.5 TeV has been applied to 
increase the event generation efficiency. For jet cluster¬ 
ing, stable particles with lifetimes greater than 3 • 10“^^ 
seconds are selected, and neutrinos are excluded. 

The SM background predictions were performed at 
leading order (LO) in QCD with the MadgraphS pro¬ 
gram. We include the following irreducible background 
processes in this study; (i) VVjj —>■ 4rjj produc¬ 
tion {V = Z,^*) via purely electroweak couplings, (ii) 
ZZjj ^Tjj production via the presence of the strong 
coupling in the Feynman amplitudes, and (iii) HHjj —>■ 
4rjj production via purely electroweak couplings. The 
VV production from the vector boson scattering topol¬ 
ogy was computed separately from ZZ events where the 
jets are radiated from QCD vertices. The interference 
between the ZZjj amplitudes with and without QCD 
vertices has been shown to be about 7% in the relevant 
phase space at the LHC [72]. We consider this inter¬ 
ference effect negligible for the purposes of this study. 
We have also neglected HHjj production via gluon fu¬ 
sion as this contribution is suppressed by the selection 
requirements favoring the VBS topology. The transverse 
momenta of the generated r leptons were required to be 
Pt{t) > 100 GeV, with the leading r lepton required to 
have pt{t) > 300 GeV, and all r leptons are required to 
have pseudo-rapidity l77(T)j < 3, to increase the efficiency 
for event generation. 

Next-to-leading order QCD /c-factors have been esti¬ 
mated [73, 74] to be ~ 50% for QCD production of VVjj 
and < 10% for VBS. As the latter background is domi¬ 
nant in this analysis, the relevant NLO-QCD correction 
is both small and similar for the signal VBS topology. 



(a)The pseudo-rapidity distributions of the forward jets. 



(b)The pseudo-rapidity distribution of the forward jet with larger 
absolute pseudo-rapidity. 

FIG. 4. The pseudo-rapidity distributions of the two for¬ 
ward jets. Generator-level cuts mjj > 1.5 TeV and PT(jet) > 
50 GeV have been applied on the jets for the samples shown. 
Furthermore, the generator-level cuts on the samples in all 
figures include pt{t) > 100 GeV, > 300 GeV and 

|77(t)| < 3. The signal distribution in all figures corresponds 
to rrin = 3 TeV with a fractional width of 20%. 
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Thus, NLO-QCD corrections are expected to have a neg¬ 
ligible impact on the results. 

A potential additional source of signal is the ZZ decay 
mode of the 77 resonance, via the 77 —>■ ZZ —^ 4r chan¬ 
nel. Due to the factor of ^ 2 smaller branching ratio 
for Z ^ TT compared to H ^ rr, this contribution is 
small compared to the signal we have considered. For 
simplicity we have neglected this contribution, yielding 
conservative results for signal sensitivity. Alternatively, 
the T] ^ ZZ ^ At channel may be distinguished from 
the rj ^ HH ^ At channel using advanced analysis 
techniques, and combined with the Z —>■ ee//r/r chan¬ 
nels to check the branching ratios of the resonance to the 
Goldstones. 


IV. KINEMATIC DISTRIBUTIONS 

Figure 4 shows the pseudo-rapidity distribution of the 
forward jets and the distribution of pseudo-rapidity of the 
forward jet with larger absolute pseudo-rapidity. For¬ 
ward jets are defined as the jet pair with the largest 
invariant mass rrijj in the event. Studies at the LHC 
have shown this criterion to be effective in identifying 
the forward tagging jets in VBS. The pseudo-rapidity 
of the more-forward jet gives an indication of the re¬ 
quired rapidity coverage of the detector. The figures 
show that processes involving QCD vertices produce jets 
which are more central, while the VBS topology for SM 
and rj resonance processes tends to produce jets in the 
forward direction. Furthermore, the 77 resonance process 
is mediated by longitudinal VBS, where the scattered 
quarks emerge at higher rapidities as compared to the 
SM VBS process mediated predominantly by transverse 
vector bosons. Thus, new physics that is specific to lon¬ 
gitudinal vector bosons will be a primary driver for max¬ 
imizing forward rapidity coverage. 

Figure 4 motivates the requirement | 77 (jet)| > 2 to 
suppress diboson production via QCD processes with¬ 
out significant loss of signal efficiency, hence we apply 
this generator-level cut for the following distributions and 
studies. QCD showering by pythia can generate addi¬ 
tional jets in the central region. The distribution of the 
difference in pseudo-rapidities Aijjj of the forward jets is 
shown in Fig. 5. This distribution motivates the addi¬ 
tional selection requirement of \^rijj\ > 5 after pythia 
showering. 

The scalar sum Ht = DIptI of all visible objects with 
pseudo-rapidity I 77 I < 2.5 in the event, called “central 
Ht" {Ht), is sensitive to QCD radiation accompanying 
the bosons and jets in the QCD-induced ZZ+2j process. 
As shown in Fig. 5, the requirement < 300 CeV 
suppresses this background. 

Figure 6 shows the distribution of the pt of the for¬ 
ward jets, and their mjj distribution. The pt spectrum 
of the forward jets emitting longitudinal vector bosons 
due to the signal process VlVl —>■ 77 is softer than the 
corresponding spectra from electroweak VBS and QCD 


processes, motivating the cut pT(jet) < 500 CeV to sup¬ 
press the latter backgrounds. 

The T leptons from the p ^ HH ^ At decay are 
produced more centrally and with higher pT than the 
backgrounds. We find that good coverage of r-jets up 
to 1 77 1 < 3 is adequate to have high acceptance for the 
signal, as shown in Fig. 7. The inclusive pt spectrum of 
all r-jets is also shown in Fig. 7. The pt distributions of 
the leading and the next-to-leading r-jets, ranked in pT, 
are shown in Fig. 8 . 

The vector sum = l^prl of all detected objects, 
which defines the missing transverse energy, is sensitive 
to the angular correlations between the r-neutrinos emit¬ 
ted in the decay of the r leptons. The r-leptons from 
Higgs boson decays have anti-parallel spins in the Higgs 
rest frame due to the Higgs being a scalar boson. Civen 
the V — A nature of the r-lepton decay vertex, the r neu¬ 
trinos are preferentially emitted parallel to each other in 
the rest frame, increasing the ^t. In comparison, the 
T neutrinos are emitted anti-parallel to each other in Z 
boson rest frame due to the unit spin of the latter, prefer¬ 
entially reducing the The distribution of is shown 
in Fig. 9. 

Additional discrimination between signal and back¬ 
ground processes is provided by the invariant mass of 
combinations of r-jets. We combine the r-jets from a 
given Higgs or Z boson decay, and average the two re¬ 
sulting invariant masses in the event. The average re¬ 
constructed boson mass distribution is shown in Fig. 10. 
The peak of the reconstructed Higgs boson mass is shifted 
to a higher value compared to the reconstructed Z bo¬ 
son mass, as expected. The use of sophisticated mass- 
reconstruction techniques that have been developed for 
di-r resonances may be used to recover information lost 
with the neutrinos, which may compensate for the exper¬ 
imental resolution on the visible momenta. Also shown 
in this figure are the distributions of the invariant mass 
all T-jets, and the combination of all r-jets and setting 
the ^z= 0. 


V. RESULTS 

After applying the selection cuts, the dominant irre¬ 
ducible background with the 4r final state is ZZ —)■ At 
production in the VBS topology. We combine the in¬ 
formation in the following distributions: the pt of the 
forward tagging jets and their pseudo-rapidity separa¬ 
tion, the Pt of the leading and sub-leading r jets, the 
^T, the iJf., the average di-r mass, the all-r mass and 
the all-T-|-^T mass, using a Boosted Decision Tree (BDT) 
algorithm to separate the 77 —>■ iJiJ —>■ 4 t signal from 
the VBS ZZ —s. At background. The resulting distri¬ 
butions of the BDT score for the signal and this domi¬ 
nant background are shown in Fig. 12. We quantify the 
discovery reach for the signal by computing the quan¬ 
tity CLh = P(Q < Qobsll^), the probability for the test- 
statistic Q to be smaller than the observed value given the 
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background-only hypothesis. When 1 — CLb < 2.8 x 10“^ 
the background-only hypothesis is rejected at 5cr signif¬ 
icance. The 5cr-discovery mass reach for the rj ^ HH 
resonance for different fractional widths and integrated 
luminosities is shown in Table I. 


TABLE I. 5cr discovery mass reach for the 77 —>■ —>■ 4r 

resonance, at a pp collider with ^/s = 100 TeV, as a function 
of integrated luminosity jC. 


c 

(ab- 1 ) 

r/M = 5 % 

TTi,, (TeV) 
r/M = 20 % 

r/M = 70% 

1 

0.85=- 

1.75 

2.81 

3 

1.33 

2.25 

3.42 

10 

1.78 

2.90 

4.18 

30 

2.30 

3.56 

4.94 

100 

2.90 

4.33 

5.83 


^ The minimum px cuts on the r lepton have been reduced for 
this mass point. 


Table II shows the dependence of the 5cr-discovery 
mass reach on the minimum pp cut applied on the for¬ 
ward tagging jets. The mass reach reduces by about 22% 
for every 20 GeV increase in the prOet) cut. Thus it is 
beneficial to maintain as low a prOet) cut as possible. 
Similarly, the forward rapidity coverage of jets is impor¬ 
tant. The dependence of the resonance mass reach as a 
function of the maximum jet rapidity detectable is shown 
in Table III. Coverage up to jet rapidity of 6-7 is desirable 
for a 100 TeV pp collider. 

TABLE 11. 5a discovery mass reach for the 77 —>■ HH —>■ 4r 
resonance, at a pp collider with ^/s = 100 TeV and C — 
10 ab~^, for various cuts values on minimum px of the forward 
jets. The fractional width of the 77 resonance is set to T /M = 
20 %. 


Pt'" (GeV) 

30 

50 

70 

90 

110 

rtir, (TeV) 

3.53 

2.90 

2.35 

1.92 

1.56 


TABLE III. 5(7 discovery mass reach for the p ^ HH —>■ 4r 
resonance, at a pp collider with yfs = 100 TeV and C = 
10 ab“^, for various cuts values on the maximum rapidity (p) 
of the forward jets. The fractional width of the p resonance 
is set to r /M = 20%. 


max 

y 

8 

7 

6 

5 

4 

rur, (TeV) 

2.9 

2.9 

2.81 

2.42 

1.75 


Tables IV and V summarize the discovery mass reach 
as a function of integrated luminosity (£) and col¬ 
lider center-of-mass energy ^/s, for resonance widths of 
20% and 70% respectively. The corresponding results are 
also presented in Figs. 13 and 14. The figures show that 
over a wide range of C and y/s, can be described 
fairly well by a power-law dependence on £, 

oc £“ (3) 


TABLE IV. 5(7 discovery mass reach for the 77 —>■ HH —>■ 
4r resonance, as a function of the y/s of a pp collider. The 
fractional resonance width r,,/ 77 i,, is fixed at 20%. These 
results are illustrated in Fig. 13. 


£ 


7 T 7 ,, (TeV) 


(ab-i) 

Vs = 50 TeV 

yi = 100 TeV 

Vi = 200 TeV 

1 

1.26 

1.75 

2.27 

3 

1.58 

2.25 

2.88 

10 

2.02 

2.90 

3.66 

30 

2.49 

3.56 

4.44 

100 

3.06 

4.33 

5.38 

TABLE V, 

. 5(7 discovery 

mass reach for the 

p -7- HH 4r 

resonance, 

as a function of the Vi of a pp collider. The frac- 

tional resonance width Tn/rriri is fixed at 70%. These results 

are illustrated in Fig. 14. 



£ 


7 T 7 ,, (TeV) 


(ab-i) 

yfs = 50 TeV 

yi = 100 TeV 

Vi = 200 TeV 

1 

1.89 

2.81 

3.85 

3 

2.31 

3.42 

4.65 

10 

2.83 

4.18 

5.63 

30 

3.36 

4.94 

6.60 

100 

3.97 

5.83 

7.74 

where the 

a values are 

independent of -^s. We find that 

the results of tables IV and V can be parameterized by 

a = 0.20 (0.16) for a 

resonance width of 20% (70%). 


Equivalently, an increase in £ by a factor of 10 raises the 
discoverable mass by 58% (45%). The gain in mass reach 
with C is slightly more rapid at low C and slightly slower 
at high £. 

We also attempt to describe the dependence of on 
\/s by a power law, 

oc (v^)^ (4) 

and find that (3 values are fairly independent of the inte¬ 
grated luminosity. If ^/s is increased from 50 TeV to 100 
TeV, /3 = 0.50 (0.56) for a fractional resonance width of 
20% (70%). If is increased from 100 TeV to 200 TeV, 
the corresponding value of the power /3 = 0.34 (0.43) fits 
the results. Thus, the scaling behavior for the discov¬ 
ery mass reach as a function of y/s can be approximated 
by a power-law behavior but the sensitivity gain starts 
to saturate at the higher collider energies. A reasonable 
approximation is obtained by fitting the gain from 50 
TeV to 200 TeV, yielding /3 = 0.42 (0.50) for T/M =20% 
(70%). Equivalently, a doubling of the collider energy 
increases the discovery mass reach by 33-40% depending 
on the resonance width, with a somewhat larger (smaller) 
increase at lower (higher) energies. 

It is interesting to evaluate the trade-off between col¬ 
lider energy and integrated luminosity for a given discov¬ 
ery mass reach. For a fractional resonance width of 20% 
(70%), a factor of two in collider energy is equivalent 
to a factor of 4.3 (8.7) in integrated luminosity. For a 



























7 


narrow (implying weakly coupled) resonance, integrated 
luminosity is more effective while for a wide (implying 
strongly coupled) resonance, collider energy is more ef¬ 
fective as a means of increasing the mass reach. 

Studies of the sensitivity of the HL-LHC for resonances 
in vector boson scattering indicate a discovery potential 
of mres/ffres of ~ 500 GeV [75], where g^es is the resonance 
coupling. While the resonance model and decay channel 
used in that study were different from ours, the sensitiv¬ 
ity from Table I is about a factor of 4-7 higher than that 
at the LHC. Comparative studies between the HL-LHC 
and a 100 TeV pp collider [76] regarding the sensitivity 
to dimension-8 operators in vector boson scattering can 
also be interpreted as a factor of four higher mass scale 
being probed at the 100 TeV collider. 


VI. CONCLUSIONS 

The sensitivity to a resonance in longitudinal VBS in 
the mass range of 1.5-5 TeV decaying to HH —>■ 4t is 
discussed for a 100 TeV pp collider. In a benchmark 
model motivated by the spontaneous breaking of a global 
SO{5) symmetry to S'0(4), the SU{2)l Higgs doublet 
field contains the set of four Goldstone modes, which are 
derivatively coupled to this resonance. The resonance 
decays democratically to longitudinal and Z bosons 
and Higgs bosons with 2:1:1 proportion for the respec¬ 
tive branching ratios. We have used leading-order cross 
sections for the signal and background processes, as the 
dominant background of VV production {V = Z,j*) in 
the VBS topology is also purely electroweak and higher- 
order QCD corrections can be expected to increase the 
signal and background cross section in comparable pro¬ 
portion. 

A reconstruction efficiency for r leptons of 60% and 
a corresponding QCD jet efficiency of a few percent is 
assumed in this study. Future detectors are expected 
to maintain this r reconstruction performance, currently 
achieved by the LHC experiments, at higher transverse 
momenta with acceptance up to a pseudo-rapidity of 3. 

This study shows that, for an integrated luminosity 
of 10 ab“^ at a pp collider with ^/s = 100 TeV, a 5cr 
discovery reach of 2.90 TeV for the mass of the reso¬ 
nance can be achieved, assuming a width of 20%. For 
widths varying between 5% and 70%, the corresponding 
mass reach varies from 1.78 TeV and 4.18 TeV. For a fac¬ 
tor of three increase in integrated luminosity, the mass 
reach increases by 25% (19%) for fractional resonance 
width of 20% (70%), and this power-law scaling behav¬ 
ior is independent of the collider energy. An approximate 
power-law scaling dependence on ^/s is also found, where 
a doubling of the collider energy increases the discovery 
mass reach by 33-40% depending on the resonance width, 
with a larger proportionate increase at lower energies. 

The trade-off between collider energy and integrated 
luminosity for a given discovery mass reach favors lumi¬ 
nosity for a narrow, weakly-coupled resonance and energy 


for a wide, strongly-coupled resonance. For a fractional 
resonance width of 20% (70%), a factor of two in collider 
energy is equivalent to a factor of 4.3 (8.7) in integrated 
luminosity. 

New physics in longitudinal vector boson scattering 
puts strong requirements on the detection of forward jets 
at high rapidities and low pT- The mass reach is reduced 
by about 22% for every 20 GeV increase in the mini¬ 
mum pT(jet) requirement on the forward tagging jets. A 
minimum rapidity coverage for the forward jets of 6-7 is 
desirable. For p^Oet) > 50 GeV, the mass reach is re¬ 
duced by 14% if the rapidity coverage is reduced from 6 
to 5. Further reducing the rapidity coverage to 4 causes 
the mass reach to drop by another 28% relative to the 
coverage of 5. If VBS jets could be distinguished from 
pileup jets at lower pT, the gains from extended rapidity 
coverage would be even higher. 

This paper also highlights the importance of the H —>■ 
rr decay channel as a relatively clean mode of identifica¬ 
tion for double-Higgs production. Distinguishing highly- 
boosted r leptons withpr(T) ^ 1 TeV from QCD jets and 
electrons presents a challenge that could be addressed 
by high-granularity electromagnetic calorimeters. Used 
in conjunction with tracking detectors having good two- 
track resolution, such calorimeters could measure the in¬ 
dividual charged particles and photons within the r-jets 
with sufficient spatial and energy resolution. Resolv¬ 
ing the substructure in r-jets and possibly extracting 
information on r polarization could allow future detec¬ 
tors to maintain or even surpass the r identification ef¬ 
ficiency and background rejection that has been demon¬ 
strated by the LHC experiments at lower transverse mo¬ 
menta. Thus, the H ^ tt mode could provide a discov¬ 
ery channel for resonant double-Higgs production with 
good signal-to-background ratio, albeit with low statis¬ 
tics, similar to the H —>■ ZZ M {i = e, p) channel for 
the Higgs boson discovery at the LHC. The fact that the 
H ^ TT branching ratio is approximately equal to the 
combined branching rato for Z —>■ ee, pp supports our 
emphasis on r detection. 
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(a)The distribution of the difference Arjjj in pseudo-rapidities of the 
two forward jets. 



(b)The distribution of the scalar sum of the pT of the central jets. 

FIG. 5. The samples shown in this figure include the 
generator-level cuts mjj > 1.5 TeV, and pT(jet) > 50 GeV. 

The distributions motivate the additional requirements of 
\^Vjj \ > 5 and i/y < 300 GeV. 


(a)The pT distribution of forward jets. 



(b)The distribution of dijet mass of the forward jets. 

FIG. 6. The Pt and dijet mass distributions of the two 
forward jets. The samples shown in this figure and subse¬ 
quent figures include the generator-level cuts rrijj >1.5 TeV, 
50 < PT(jet) < 500 GeV, |7j(jet)| > 2, \Ar]jj\ > 5 and 
H!r < 300 GeV. 
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FIG. 



(a)The pseudo-rapidity distribution of r-jets. 



(b)The pt distribution of r-jets. 

7. The pseudo-rapidity and pt distributions of r-jets. 



(a)The pT distribution of the leading r-jet. 



(b)The Pt distribution of the next-to-leading r-jet. 


FIG. 8. The pr distributions of the highest-pT r-jets. 
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FIG. 9. The distributions of 



FIG. 10. The distribution of the per-event average of the 
reconstructed masses of the r-jet pairs from V ^ tt {H —>■ 
rr) in VVjj {HHjj) events, where V = 7*/^. 



(a)The distribution of invariant mass of all reconstructed r-jets. 



(b)The distribution of invariant mass of all reconstructed r-jets and 
^T, setting F/z = 0. 

FIG. 11. The invariant-mass distributions involving multiple 
reconstructed r-jets. 
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FIG. 14. 5cr discovery mass reach for the ry —>■ HH —>■ 4r 
resonance, as a function of the integrated luminosity and y/s 
of a pp collider. The fractional resonance width is 

fixed at 70%. 


FIG. 12. The distributions of the BDT score for the rj —> 
HH —>■ At signal with m,, = 3 TeV, and the VBS VV —>■ 4t 
background, where V = jZ. 



FIG. 13. 5cr discovery mass reach for the p ^ HH ^ At 
resonance, as a function of the integrated luminosity and ^/s 
of a pp collider. The fractional resonance width Vnlrrin is 
fixed at 20%. 
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Decay channel 


Branching ratio 


Uncertainty 


bbbb 

rrbb 

W+i-^ lv)W-{-^ lv)bb 

TTTT 
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W+{^ lu)W-{^ lu)Tr 

77rr 

bbfi^ fj,~ 

z\-^ l+l-)Z{^ l+l-)bb 
bbZ{-^ i+r)y 

W+i-^ lv)W-{-^ lu)W+{^ lu)W-{^ Iv) 
77 lU+(-> lv)W-{-^ Iv) 

TTfJ,^ fj,~ 

Z{-^ l+l-)Z{^ l+l-)Tr 
7777 

W+i-^ ly)W-{-^ lv)li+tj.- 

Z{-^ l+l-)Z{^ l+l-)W+{^ lv)W-{-^ Iv) 

W+{-^ lv)W-{-^ lv)Z{^ l+l-)y 

77/i+^" 

yyZ{^ l+l-)Z{-^ 1+1-) 
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Z{-^ l+l-)jfi+fi- 

z{^ i+i-)z{^ i+i-)z{-^ i+rh 
z\-^ i+i-)z{^ i+i-)z{^ i+i-)z{^ i+r) 

Z{^ l+r ^-yZi^ l+l-)'y 


3.33 ■ 

10-1 

±2.20- 

10-^ 

7.29 • 

10"^ 

± 4.80 • 

10"® 

1.09 ■ 

10-^ 

± 5.93 • 

10-1 

3.99 • 

10"® 

±4.55 • 

10-1 

2.63 • 

10-^ 

± 1.58 • 

10-1 

1.20 • 
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4.78 • 

10-" 

± 4.92 • 

10-® 

5.35 • 

10"® 

± 3.95 • 
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4.59 • 

10"® 
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10-® 
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10"® 
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10-® 

1.49 • 

10"® 
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TABLE VI. Branching ratios for final states arising from double-Higgs prodnction, with the reqnirement of leptonic decays of 
W and Z bosons. 
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FIG. 15. An event display showing the 77 resonance (m^ = 3 TeV) produced via longitudinal vector boson scattering, and 
decaying to 1/11 —>■ 4r at a 100 TeV pp collider, (left) In the transverse view, two forward jets with transverse momenta above 
50 GeV are shown with the small yellow cones, and two jets arising from H —>■ tt are shown with the large yellow cones, (right) 
A 3D view of the same event. The jets are reconstructed using the anti-fer algorithm [22] using the FastJet package [23]. The 
event display was created using the Delphes fast simulation [56], HepSim [59] and the Snowmass detector setup [57]. The blue 
lines show charged hadrons. See the text for details. 





